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ABSTRACT 

We report on an observation of the neutron star low-mass X-ray binary Serpens X-l, made with NuSTAR. The 
extraordinary sensitivity afforded by NuSTAR facilitated the detection of a clear, robust, relativistic Fe K emission 
line from the inner disk. A relativistic profile is required over a single Gaussian line from any charge state of 
Fe at the 5cr level of confidence, and any two Gaussians of equal width at the same confidence. The Compton 
back-scattering “hump” peaking in the 10-20 keV band is detected for the first time in a neutron star X-ray binary. 
Fits with relativistically blurred disk reflection models suggest that the disk likely extends close to the innermost 
stable circular orbit (ISCO) or stellar surface. The best-fit blurred reflection models constrain the gravitational 
redshift from the stellar surface to be jns ^ 0.16. The data are broadly compatible with the disk extending to the 
ISCO; in that case, zns ^ 0.22 and Rns ^ 12.6 km (assuming Mns = 1.4 M 0 and a = 0, where a = cJ /GM 2 ). If 
the star is as large or larger than its ISCO, or if the effective reflecting disk leaks across the ISCO to the surface, the 
redshift constraints become measurements. We discuss our results in the context of efforts to measure fundamental 
properties of neutron stars, and models for accretion onto compact objects. 
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1. INTRODUCTION 

The equation of state of ultradense matter remains poorly 
known (see, e.g., Lattimer 2011). Such conditions cannot be 
replicated in terrestrial laboratories, leaving only astrophysical 

studies of neutron stars to solve this fundamental problem. The 

mass-radius relation for neutron stars depends sensitively on the 

equation of state of the ultradense stellar interior. Neutron star 
radius measurements are a particularly powerful discriminant, 

although extreme mass measurements can also help (e.g., 
Demorest et al. 2010). X-ray observations hold tremendous 
promise for radius measurement. Current instruments can easily 

measure thermal spectra from X-ray bursters (e.g., Guver et al. 

2010) and quiescent X-ray transients (e.g., Guillot et al. 2013), 

although systematic uncertainties involving source distance, 

atmospheric modeling, and additional spectral components can 
complicate the extraction of reliable stellar radii (Boutloukos 

et al. 2010; Guver et al. 2012a, 2012b; Steiner et al. 2013). 

Whenever an accretion disk is illuminated by an external 

source of hard X-rays, disk reflection is expected. The emergent 
reflection spectrum is shaped by Doppler shifts and gravitational 
redshifts; the extremity of these shifts can be used to measure 
the location of the inner disk with respect to the compact object 
(for a review, see Miller 2007; also see Miller et al. 2009). 
The most prominent parts of a disk reflection spectrum are the 


Fe K emission line and the Compton back-scattering “hump” 
that peaks between 10-30 keV. In principle, disk reflection 
spectroscopy has the potential to constrain the radius of a 
neutron star, since the disk must truncate at the surface, if not at 
larger radii (e.g., Cackett et al. 2008, 2010). 

Spectra obtained with XMM-Newton and Suzaku have re- 
vealed potentially relativistic line shapes and have started to 
illustrate the power of disk reflection spectra to constrain the 
fundamental parameters of neutron stars (e.g., Bhattacharyya 
& Strohmayer 2007; Cackett et al. 2008, 2010; Di Salvo et al. 
2009; Iaria et al. 2009; D’Ai et al. 2009; Papitto et al. 2009; 
Miller et al. 2011; Egron et al. 2013). A difficulty with X-ray 
CCD detectors is that pile-up can distort the shape of Fe K lines. 
Indeed, Ng et al. (2010) suggested that Fe K lines in neutron 
stars may actually be narrow and symmetric, and merely appear 
to be relativistic owing to pile-up distortions. 

However, this possibility is not borne out by simulations (see, 
e.g., Miller et al . 20 1 0) , which instead suggest that pile-up acts to 
falsely narrow lines. Data obtained with gas spectrometers (e.g., 
Lin et al. 2010; Cackett et al. 2012), which do not suffer from 
photon pile-up, also point to broad lines, but such studies are 
clouded by low spectral resolution. Thus, uncertainties remain. 
The detection of a Compton back-scattering “hump” would also 
help to resolve the origin of Fe K lines in neutron stars, but no 
detection has previously been achieved. 
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Motivated by the ability of NuSTAR (Harrison et al. 2013) 
to reveal Fe K lines and broad-band disk reflection spectra in 
black holes (e.g., NGC 1365, GRS 1915+105, and Cygnus 
X-l; Risaliti etal. 2013; Miller etal. 2013; Tomsicket al. 2013), 
free of photon pile-up distortions, we made a pilot observation 
of the well-known neutron star X-ray binary and “atoll” source 
Serpens X-l. 

2. OBSERVATIONS AND DATA REDUCTION 

NuSTAR observed Serpens X-l twice on 2013 July 12 
and 13 (Sequence IDs 30001013002 and 30001013004). The 
data were screened and processed using the NuSTAR Data 
Analysis Software version 1.1.1, resulting in a total of 43.4 ks 
of on-target time. The source is relatively bright, so the net 
exposure after instrumental deadtime is taken into account is 
30.4 ks. Spectra from the FPMA and FPMB detectors were 
extracted from 120" regions centered on the source position, 
using the “nuproducts” FTOOL. Corresponding response files 
appropriate for the pointing type (on-axis), source type (point), 
and extraction regions were also created. Backgrounds were 
created using the “nuskybgd” tool, which self-consistently 
accounts for background variations across the face of the 
detector. 

The spectra were analyzed using XSPEC version 12.8 
(Arnaud & Dorman 2000). The y 2 statistic was used to mea- 
sure the relative quality of different spectral models. Following 
Miller et al. (2013), “Churazov” weighting was employed in all 
fits to govern the influence of bins with significant but progres- 
sively less signal at high energy (Churazov et al. 1996). Errors 
are reported as the la confidence level on all measurements. 
All fits were made assuming a neutral equivalent neutral hy- 
drogen column density of 4.0 x 1 0 21 cm 2 (based on Dickey & 
Lockman 1990), via the “tbabs” model, setting “vern” cross sec- 
tions (Verner et al. 1996) and “wilm” abundances (Wilms et al. 
2000). 

3. ANALYSIS AND RESULTS 

Serpens X-l is known to exhibit Type-1 X-ray bursts at a 
rate of approximately 0. 1 per hour (Galloway et al. 2008). Light 
curves were therefore created via the “nuproducts” tool and 
were searched for X-ray bursts; no bursts were detected. We 
also constructed an X-ray color-color diagram using energy 
selections typical of diagrams made using RXTE data, in order 
to provide the best correspondence. The NuSTAR observation 
was found to sample the usual “banana” branch rather than the 
“island” state or any other period of extreme or unusual behavior 
(see van der Klis 2006 for a description of atoll sources). 

An initial inspection and comparison of the FPMA and FPMB 
spectra suggested that the source is confidently detected out to 
40 keV. We therefore restricted subsequent spectral analysis to 
the 3 — 40 keV band. Initial fits revealed a striking agreement 
between the FPMA and FPMB spectra: when a constant was 
floated between the spectra to account for uncertainties in the 
flux calibration of the detectors, a value of 1.001 was mea- 
sured. We therefore created a single combined total spectrum, 
background spectrum, and ancillary response matrix using the 
FTOOL “addascaspec.” A single redistribution matrix file was 
created via “addrmf,” weighting the single responses by their 
relative exposure times. The analysis described below pertains 
to the single combined spectrum. 

A model consisting of a disk blackbody component 
(“diskbb,” see Mitsuda et al. 1984), a simple single-temperature 
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blackbody (“bbody”), and a power-law provides a good de- 
scription of the time-averaged spectrum. This model can be 
interpreted fairly simply in terms of emission from the accre- 
tion disk, the stellar surface, and non-thermal emission that 
may arise through Comptonization and/or magnetic processes. 
The simple blackbody may not represent direct emission from 
the stellar surface, nor from the entire surface. It may rep- 
resent emission that has been modified by a boundary layer 
that can also be modeled via low-temperature, optically thick 
Comptonization. However, that prescription is similar to a hot 
blackbody. Although this three-component model describes the 
continuum shape fairly well, it is not a formally acceptable fit: 
X 2 /v = 3282.8/919 = 3.57 (also see Figures 1 and 2). 

The fit is greatly improved when a relativistic line function 
is added: / 2 /u = 1010.9/912 = 1.11 (see Figure 1). This 
model fits the line with the “kerrdisk” function (Brenneman 
& Reynolds 2006). This line model was employed because it 
allows for small changes to the spacetime induced by spin. The 
innermost stable circular orbit (ISCO) radius changes slowly for 
low values a, and neutron stars in X-ray binaries have low spin 
parameters (a <( 0.3, where a = cJ/GM 1 ; see, e.g., Galloway 
etal. 2008; Miller etal. 2011). 

We measure a disk blackbody temperature of kT = 
1.833(5) keV and a normalization of K = 29.0(3), a 
simple blackbody temperature and normalization of kT = 
2.495(5) keV and K = 0.0222(3), and a power-law 
index and normalization of T = 3.41(3) and K = 

1 ,98(4)photons cm 2 s 1 at 1 keV. For “kerrdisk,” we fixed outer 
emissivity index to equal the inner emissivity index, constrained 
the energy to the range possible for Fe (6.40-6.97 keV), and 
measured the line centroid energy, spin parameter, inclination, 
inner emission radius, and line flux. We measure a line cen- 
troid energy of E = 6.97_o.oi keV (consistent with H-like Fe 
XXVI), an emissivity index of q = 3.5(1), a spin parameter 
of a — 0.17(2), an inner disk inclination of i = 18(2)°, an 
inner radius of r m = 1.8(1)/tsco (10.6 ± 0.6 r g ), and a flux 
of 5.5 x 10~ 3 photonscm -2 s -1 . This flux translates into a line 
equivalent width of W = 91(2) e V. 

Fits with a relativistic line function are strongly required over 
a Gaussian model for the emission line. The best fits with the 
“kerrdisk” and “diskline” functions give y 2 /v — 1010.9/912. 
Instead modeling the line with a Gaussian corresponding to any 
charge state of Fe (6.40-6.97 keV) gives y 1 /v — 1052.1/915 
( E — 6.53(2) keV, a = 0.44(2) keV). Via an F-test, a relativistic 
line function is therefore required at the 5.3a level of confidence. 

It is also possible to test whether or not the line profile 
is actually a composition of different plausible charge states. 
Adding two Gaussians from any charge state of Fe, constrained 
to have the same width but allowed to have different strengths, 
gives a significantly worse fit than a single relativistic line: 
y 2 /v — 1266.1/915. This effectively eliminates the possibility 
that the line profile arises in any form of diffuse gas far from the 
stellar surface. 

Using this phenomenological model, we measure an unab- 
sorbed flux of 6. 1 x 10 _9 erg cm~ 2 s -1 in the 3.0-40.0 keV band, 
implying 1 .5 x 10~ 8 erg cm~ 2 s 1 in the 0.5-40.0 keV band. The 
3-40 keV range is closer to that common in X-ray monitors 
and corresponds to 0.2-0. 3 Crab. The distance to Serpens X-l 
is likely 7.7 ± 0.9 kpc (based on X-ray burst properties; see 
Galloway et al. 2008). A luminosity of 1.1(2) x 10 38 ergs _I 
follows from the inferred 0.5-40.0 keV flux. It is likely that 
Serpens X-l was accreting at a high fraction of its Eddington 
limit during our observation. 
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Figure 1. NuSTAR 3 — 40 keV spectrum of Serpens X-l is shown above, with the FPMA and FPMB detectors combined. Left: the best-fit phenomenological spectral 
model is shown here. The model consists of disk blackbody (blue), simple blackbody (red), power-law (cyan), and relativistic disk line (black) components. The edge 
feature between 9-10 keV is consistent with the K edge of Fe XXVI, and the flux excess in the 10-20 keV band is consistent with the Compton back-scattering 
“hump” expected from disk reflection. Right: the best relativistically blurred disk reflection fit (Model 5 in Table 1) is shown here. The line function has been replaced 
with a version of “reflionx” suited to blackbody illumination of the disk, and convolved with the “kerrconv” function. A significantly better fit is achieved. 

(A color version of this figure is available in the online journal.) 



Figure 2. Relativistic Fe K disk line in the combined FPMA and FPMB spectrum 
of Serpens X- 1 is shown above. The best three-component continuum model was 
fit over the 3-40 keV band, ignoring the 4. 5-7. 5 keV range to avoid biasing the 
continuum fit. The best “kerrdisk” model, derived from fitting the full bandpass 
including the Fe K region, is shown in red. The spectrum was binned for visual 
clarity only. 

(A color version of this figure is available in the online journal.) 

Figure 1 shows clear hallmarks of disk reflection, even after 
fitting the relativistic line, including an edge between 9-10 keV 
consistent with H-like Fe XXVI, and a flux excess in the 
10-20 keV range consistent with a Compton back-scattering 
hump. It is also clear in Figure 1 that most of the flux capable 
of ionizing Fe comes from the blackbody component. To 
self-consistently model the reflection, then, we employed a new 
version of the well-known “reflionx” model (Ross & Fabian 
2005), wherein a blackbody illuminated the disk instead of the 
power-law. 

The “kerrdisk” line model was replaced by “reflionx,” con- 
volved with the kernel of “kerrdisk,” called “kerrconv.” Its pa- 
rameters include inner and outer disk emissivity indices (c/ OLU 
was fixed to q m in our fits), the inclination of the inner disk, the 


spin parameter of the compact object, and the inner and outer 
emitting radii (in units of the ISCO radius; R out — 400 was fixed 
in all cases). The “reflionx” model variables include the ioniza- 
tion parameter of the disk (£ = L/ nr 2 ), the abundance of iron 
( Af e ), the blackbody temperature (tied to that of the independent 
blackbody), and the flux. 

Table 1 lists the results of models for the time-averaged 
spectrum of Serpens X- 1 , including relativistically blurred disk 
reflection (Models 1-18). The spin of Serpens X-l is unknown 
as kHz quasi-periodic oscillations and burst oscillations have not 
been detected. For completeness we considered fits with three 
spin values: a = 0.0, 0.12, 0.24. The latter two values span the 
range of spin parameters implied in other “atoll” sources (see, 
e.g., Galloway et al. 2008), assuming R\% = 10 km, Mns = 
1 .4 M e , and / NS = (2/5 )M R 2 as per a simple sphere. Similarly, 
the iron abundance of Serpens X-l is unknown, so we explored 
fits with three values (relative to solar): Ap e = 1.0, 2.0, 3.0. 
Last, in order to test the sensitivity of the spectrum to the inner 
radius of the disk, we considered fits with variable inner radii, 
and fits with r; n = 1.0/'isco- 

The best blurred reflection fits (Models 5, 11, 17) uniformly 
give significantly better fits than models without reflection: an 
/’’-test gives an F = 96.3, or a false-alarm probability of just 
1 . 1 x 10~ 21 . Combined with the ability of the reflection models 
to naturally model the 10-20 keV flux excess in Figure 1, this 
signals that the full broad-band disk reflection spectrum has 
been detected in Serpens X-l. 

All of the models in Table 1 prefer very low inclinations; 
this corresponds well with recent optical observations that 
independently point an inclination of 10° or less (Cornelisse 
et al. 2013). The results in Table 1 also clearly signal that the 
fits are largely insensitive to the spin parameter of the neutron 
star. This is consistent with expectations: even for a = 0.24, the 
ISCO radius only changes marginally to nsco = 5.89r„ (where 
r g = GM/c 2 ). 

The data show a slight preference for an enhanced iron 
abundance in Serpens X-l. Given that the binary orbital period 
of Serpens X-l appears to be just 2 hr (Cornelisse et al. 2013), it 
is possible that the companion is evolved or degenerate, skewing 
e.g., Fe/O abundance ratios. Enhanced abundances are partly 
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Table 1 

Relativistically Blurred Disk Reflection Models 


Spectrum 

Model 

& Tdisk 

^disk 

&7bb 

Kbb 

r 

^pow 

a 

9 

r in 

i 


Ap e 

K-refl 

X 2 /v 



(keV) 


(keV) 

(10“ 2 ) 





(ISCO) 

(deg.) 

(erg cm s -1 ) 


(10” n ) 


Time-avg 

1 

1.60(1) 

45(2) 

2.33(1) 

2.52(8) 

3.17(3) 

1.3(1) 

0.0* 

2.21(4) 

1.2(1) 

oo 

370(40) 

1.0* 

4.5(1) 

995.1/914 

Time-avg 

2 

1.61(3) 

44(3) 

2.33(2) 

2.5(1) 

3.18(5) 

1.4(2) 

0.0* 

2.19(3) 

1.0* 

o +2 

—8 

330(80) 

1.0* 

5(2) 

998.6/915 

Time-avg 

3 

1.61(2) 

41(2) 

2.33(1) 

2.63(6) 

3.21(3) 

1.5(1) 

0.0* 

2.26(4) 

1.3(1) 

o +2 

—8 

220(60) 

2.0* 

5(2) 

962.3/913 

Time-avg 

4 

1.64(3) 

41(3) 

2.33(2) 

2.6(1) 

3.22(4) 

1.5(1) 

0.0* 

2.23(3) 

1.0* 

q +2 

°-8 

300(90) 

2.0* 

6(3) 

966.0/915 

Time-avg 

5 

1.64(2) 

41(3) 

2.33(2) 

2.65(5) 

3.23(4) 

1.5(1) 

0.0* 

2.29(3) 

1.3(1) 

q +2 

0-8 

200(20) 

3.0* 

7(2) 

950.9/914 

Time-avg 

6 

1.64(2) 

41(3) 

2.33(2) 

2.70(5) 

3.24(4) 

1.6(1) 

0.0* 

2.27(3) 

1.0* 

8-8 

270(90) 

3.0* 

6(2) 

955.9/915 

Time-avg 

7 

1.63(3) 

42(2) 

2.34(1) 

2.37(8) 

3.21(2) 

1.48(3) 

0.12* 

2.24(4) 

1.25(3) 

8/ 2 8 

240(30) 

1.0* 

6.0(9) 

998.4/914 

Time-avg 

8 

1.61(3) 

43(3) 

2.33(2) 

2.5(1) 

3.21(3) 

1.5(1) 

0.12* 

2.20(4) 

1.0* 

0+2 

0-8 

340(90) 

1.0* 

5(2) 

998.7/915 

Time-avg 

9 

1.63(2) 

41(2) 

2.33(1) 

2.6(1) 

3.21(3) 

1.5(1) 

0.12* 

2.24(4) 

1.3(2) 

8-8 

250(60) 

2.0* 

7(3) 

962.4/914 

Time-avg 

10 

1.64(3) 

41(3) 

2.33(2) 

2.61(9) 

3.23(4) 

1.6(2) 

0.12* 

2.22(3) 

1.0* 

8 + - 2 8 

300(90) 

2.0* 

6(2) 

966.6/915 

Time-avg 

11 

1.64(2) 

40(1) 

2.33(1) 

2.69(3) 

3.22(2) 

1.51(5) 

0.12* 

2.27(3) 

1.35(3) 

8 + - 2 8 

210(20) 

3.0* 

6.1(6) 

950.9/914 

Time-avg 

12 

1.64(3) 

40(3) 

2.33(2) 

2.69(9) 

3.24(4) 

1.6(2) 

0.12* 

2.26(3) 

1.0* 

8/| 

280(90) 

3.0* 

5(3) 

957.1/915 

Time-avg 

13 

1.64(1) 

42(2) 

2.35(1) 

2.5(2) 

3.20(2) 

1.5(1) 

0.14* 

2.22(3) 

1.4(1) 

8/ 2 s 

290(40) 

1.0* 

5.2(4) 

995.8/914 

Time-avg 

14 

1.61(3) 

44(3) 

2.34(2) 

2.5(2) 

3.19(4) 

1.5(2) 

0.14* 

2.19(3) 

1.0* 

8Z 4 8 

330(90) 

1.0* 

5(2) 

1000.0/915 

Time-avg 

15 

1.63(1) 

41(2) 

2.33(1) 

2.6(1) 

3.22(2) 

1.47(6) 

0.14* 

2.25(3) 

1.4(1) 

8/ 2 8 

220(60) 

2.0* 

7(2) 

962.2/914 

Time-avg 

16 

1.64(3) 

41(3) 

2.33(2) 

2.7(2) 

3.24(4) 

1.6(2) 

0.14* 

2.23(3) 

1.0* 

q +2 

0-8 

300(90) 

2.0* 

6(3) 

967.6/915 

Time-avg 

17 

1.65(1) 

41(1) 

2.31(3) 

2.69(4) 

3.20(2) 

1.51(4) 

0.14* 

2.27(3) 

1.41(4) 

q +2 

0-8 

210(30) 

3.0* 

6(1) 

951.0/914 

Time-avg 

18 

1.65(2) 

41(3) 

2.33(1) 

2.7(1) 

3.25(4) 

1.6(2) 

0.14* 

2.27(3) 

1.0* 

q +8 

0-8 

270(90) 

3.0* 

5(2) 

958.1/915 

High-flux 

19 

1.75(5) 

37(3) 

2.40(2) 

2.6(2) 

3.30(5) 

1.7(3) 

0.00* 

2.21(7) 

i 9+0.4 
i - z -0.2 

q+3 

0-8 

380(90) 

1.0* 

5(2) 

936.7/914 

High-flux 

20 

1.73(5) 

36(4) 

2.49(4) 

2.6(2) 

3.31(7) 

1.8(3) 

0.00* 

2.21(6) 

1.0* 

q+4 

0-8 

360(80) 

1.0* 

5(2) 

937.1/915 

High-flux 

21 

1.73(6) 

36(3) 

2.39(3) 

2.6(3) 

3.35(9) 

1.8(3) 

0.00* 

2.26(6) 

1 1 +0 - 4 

q+3 

0-8 

310(90) 

2.0* 

6(3) 

928.8/914 

High-flux 

22 

1.75(6) 

35(5) 

2.38(3) 

2.6(3) 

3.33(7) 

1.8(3) 

0.00* 

2.25(6) 

1.0* 

q+4 

0-8 

330(90) 

2.0* 

6(3) 

929.6/915 

High-flux 

23 

1.73(4) 

35(4) 

2.38(2) 

2.7(2) 

3.33(7) 

1.8(3) 

0.00* 

2.31(7) 

i 2 +0 - 3 
i ' Z -0.2 

q+4 

0-8 

330(90) 

3.0* 

5(2) 

926.4/914 

High-flux 

24 

1.76(6) 

34(4) 

2.39(4) 

2.7(3) 

3.35(8) 

1.9(4) 

0.00* 

2.29(5) 

1.0* 

q+5 

0—8 

340(90) 

3.0* 

6(3 

927.5/915 


Notes. The results of fits with relativistically blurred disk reflection components (“reflionx” convolved with ‘'kerrconv”) are presented here. The disk emissivity is 
measured by the parameter q (Jr~ q ). Parameters fixed at a given value are marked with an asterisk. The ionization of the disk is measured by the ionization parameter 
£. All fits were made assuming an equivalent neutral hydrogen column density of Nu = 4.0 x 10 21 cm 2 as per Dickey & Lockman (1990), via the “tbabs” model. 
All errors quoted in the table are ltr uncertainties. 


degenerate with the normalization of the reflection component 
and its ionization. Excellent fits are achieved regardless of the 
abundance, however, and parameters of greater interest show no 
trend with abundance. 

The models wherein the inner disk radius is allowed to float 
nominally prefer a radius that is slightly larger than the ISCO; 
however, the improvement in the fit statistic is not significant, 
and 90% and/or 3 a confidence errors include the ISCO radius. 
Fits to the time-averaged spectra may be characterized as being 
broadly consistent with a disk extending to the ISCO. 

We also examined the broad flux trends within the observa- 
tion, and identified a period of high flux for consideration (see 
Figure 3). Total FPMA and FPMB spectra, background spec- 
tra, and responses were created for “GTI2” in exactly the same 
manner as for the time-averaged spectrum. The results of fits to 
this segment of the observation are also given in Table 1 (see 
Models 19-24). 

Figure 4 shows plots of A/ 2 versus r m for models 5 and 
23, in order to illustrate the sensitivity of the spectra to the 
inner extent of the disk. Fits to the time-averaged (Model 5) 
and high-flux (Model 23) spectra point to a disk that hovers 
close to the ISCO. Model 5 prefers rj n > risco > but Model 23 
suggests that r m = risen is fully compatible with the data. The 
disk continuum is less incisive than the reflection spectrum, but 
it can serve as a broad consistency check. The disk continua 
in Table 1 imply inner radii of 13-15 km (1.03-1.20 risco or 


6. 2-7. 2 r g for Mns = 1.4M 0 and a = 0) for a canonical 
spectral hardening factor of 1.7 (Shimura & Takahara 1995). 

For a = 0, the gravitational redshift implied by the fits 
is given by z + 1 = I / v /l — 2G M / rc 2 . The best models 
for the time-averaged spectra in Table 1 (e.g.. Model 5) give 
r = 1.3risco = 7.8 r g . The neutron star must be smaller than 
the inner disk radius, giving a limit of zns S? 0.16. Within the 
high flux interval, even smaller inner disk radii are preferred, 
and errors include r = risco- If the disk extends to the ISCO, 
Zns 0.22 and ANs ^ 12.6 km assuming Mns = 1.4Mq. 
The masses of neutron stars in low-mass X-ray binaries are not 
known precisely, but 1 .4 M 0 is a reasonable average (Lattimer 
2011). If the neutron star surface itself truncates the disk (and if 
any boundary layer is small in radial extent), or if the effective 
reflection radius extends within the ISCO, then these limits 
could be regarded as measurements. 

4. DISCUSSION 

We have analyzed an initial observation of the neutron star 
low-mass X-ray binary Serpens X-l. Owing to the abilities 
of NuSTAR, an Fe K line in the spectrum is revealed to be 
skewed at a high level of statistical confidence, and the Compton 
back-scattering “hump” expected in disk reflection scenarios is 
detected in a neutron star X-ray binary for the first time. Both 
findings are robust since the FPMA and FPMB detectors do 
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Figure 3. NuSTAR light curve of Serpens X-l. The interval defined as “high flux” was selected for independent spectral analysis for comparison (see Table 1). 

(A color version of this figure is available in the online journal.) 



r in / r iSC0 


Figure 4. Figure above displays the change in the fitting statistic. A/ 2 , vs. the 
inner disk radius measured via two blurred reflection fits. Model 5, which is a 
fit to the time-averages spectrum of Serpens X-l, is shown in black. Model 23, 
which is a fit to the high-flux phase of the observation (see Section 3), is shown 
in red. Both models, like others detailed in Table 1 , strongly prefer a disk that 
is close to a likely ISCO, and are statistically consistent with the disk extending 
to the ISCO itself. For plotting clarity, a constant of 1 .0 was added to all A/ 2 
values. 

(A color version of this figure is available in the online journal.) 


not suffer from flux distortions due to effects such as photon 
pile-up. Clearly, broad Fe K lines in at least some neutron star 
X-ray binaries arise in the inner disk, and can be used to probe 
the stellar properties and inner accretion flow. 

In contrast to black holes, an ISCO is not necessarily a feature 
of accretion onto neutron stars. The stellar radius could exceed 
its own ISCO, and/or the stellar magnetic field can truncate 
the disk at much larger radii. The latter is certainly the case in 
typical X-ray pulsars, but also in transient sources such as SAX 
J1808.4— 3658 and IGR J17480-2446 in Terzan 5 (Miller et al. 
20 1 1 ). Fits to the robust data obtained with NuSTAR indicate that 
the disk likely extends to the ISCO, or hovers just above such 
a radius. The high luminosity implied during this observation 
makes it more likely that the disk is not impeded by a magnetic 
field; this is also indicated by fits to a high-flux interval of the 
observation. 


If the neutron star radius in Serpens X-l is small, this 
observation represents the strongest observational indication 
of an ISCO in a neutron star low-mass X-ray binary. In this 
case, our fits require z N s ^ 0.16 (best fit), or z ^ 0.22 and 
Rt\s ^ 12.6 km if r = /'isco, assuming Mns = 1.4M 0 . If the 
star in Serpens X- 1 is larger than its ISCO, the disk may then be 
truncated by a the stellar surface and/or a thin boundary layer. 
It is also possible that material in the inner disk may leak across 
the ISCO and be reflective between the ISCO and stellar surface. 
In either of these two scenarios, the limits placed by our fits may 
be characterized as measurements. Some equations of state may 
then be incompatible with our observations (see Lattimer 201 1). 

The statistically significant skewing of the line profile clearly 
signals that dynamics and gravitation close to the neutron 
star drive our results. Flowever, the results are not entirely 
model-independent. A disk reflection spectrum is required by 
the data. Although we achieved excellent fits with “reflionx,” 
different and/or future reflection and relativistic convolution 
models may yield slightly different results. 
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